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Abstract. Dynamical urn models, such as the Ehrenfest model, have played 
an important role in the early days of statistical mechanics. Dynamical many- 
urn models generalize the former models in two respects: the number of urns is 
macroscopic, and thermal effects are included. These many-urn models are exactly 
solvable in the mean-field geometry. They allow analytical investigations of the 
characteristic features of nonequilibrium dynamics referred to as aging, including 
the scaling of correlation and response functions in the two-time plane and the 
violation of the fluctuation-dissipation theorem. This review paper contains a 
general presentation of these models, as well as a more detailed description of two 
dynamical urn models, the backgammon model and the zeta urn model. 



Prologue 

Urns containing balls, just as dice or playing cards, are ubiquitous in writings on 
probability theory, reminding us that this branch of mathematics owes its early 
developments to practical questions arising in playing games. Dynamical urn models, 
such as the Ehrenfest model, have played an important role in the elucidation of 
conceptual problems in statistical mechanics. More recently, dynamical many-urn 
models have been investigated in the mean-field geometry. These exactly solvable 
models exhibit characteristic features of nonequilibrium dynamics referred to as aging, 
such as the scaling of the correlation and response functions in the two-time plane and 
the violation of the fluctuation-dissipation theorem. 

This paper contains a didactic introduction to urn models (Section 1), a presenta- 
tion of static and dynamical properties of many- urn models (Section 2), a reminder of 
the main characteristic features of aging (Section 3), and an overview of recent results 
on two dynamical urn models, namely the backgammon model (Section 4) and the 
zeta urn model (Section 5). 

1. Urn models 

1.1. The Ehrenfest urn model 

This model was devised by P and T Ehrenfest in 1907 in their attempt to critically 
review Boltzmann's iJ— theorem. It is defined as follows: N balls are distributed in 
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two urns (or boxes). At random times, given by a Poisson process with unit rate, a 
ball is chosen at random, and moved from the box in which it is to the other box. 

Let Ni{t) (resp. N2{t) = N — Ni{t)) be the number of balls in box number 1 
(resp. number 2) at time t. For any initial configuration, the system relaxes to 
equilibrium for infinitely long times. The equilibrium state is characterized by a 
binomial distribution of the number iVi: 



/fe,cq=7'eq(iVl=fc)=2-^^(^^j (fc = 0,...,iV), (1.1) 

reminiscent of the Maxwell-Boltzmann statistics for indistinguishable particles. 

As noticed first by Kohlrausch and Schrodinger |^ , A^i (t) can be viewed as the co- 
ordinate of a fictitious walker. The temporal evolution of the occupation probabilities 
fk{t) = V{Ni{t) = k) is governed by the master equation 

dMt)^k + l N + l-k 

-— ^/fc+i(t)H — }k-i[t) - Jk{t). (1.2) 

Indeed, a move of a ball from box number 1 to box number 2 (resp. from box number 2 
to box number 1) yields fc fc — 1 (resp. fc — > fc + 1), and occurs with a rate k/N 



(resp. {N — k)/N) per unit time. The master equation (1.2) describes a non-uniformly 
biased random walk over the integers 0, . . . , l§ Hi- The first (resp. the second) term 
in the right-hand side is absent for fc = A^ (resp. fc = 0). 

The spec trum of relaxation rates of the model can be derived by looking for 
solutions to ( |l.2| ) of the form fk{t) = (t>ke~'^*- One thus obtains Am = 2m/N, 
with m — 0, . . . , N. The static solution (to — 0) is nothing but the equilibrium 
distribution ( [l.l[ ). The occupation probabilities have an exponential convergence to 
their equilibrium values, with a relaxation time 

Tcq = 1/Al = N/2. 

There are, however, much larger time scales in the Ehrenfest model. Consider 
indeed the length of time it takes for box number 1 to get empty, i.e., the first time to 
such that Ni{to) = 0. This time depends on the initial state and on the whole history 
of the system. Its mean value Tq = (to) can be evaluated as follows. The equilibrium 
probability for box number 1 to be empty, /o.oq — 2~^, is exponentially small in the 
number of balls. It is therefore expected that the typical time needed to reach this 
very rare event scales as 

To ~ l//o,cq - 2^. (1.3) 

This result can be derived in a more rigorous way. Related matters are discussed in 
References [0, p[ . The process of emptying one of the boxes is therefore characterized 



by an exponentially large time Tq. Equation (1.3) can be recast as an Arrhenius- 
like law: 

To^e^°, (1.4) 

where 5*0 = A^ In 2 is the entropy difference between the equilibrium state of the 
model, where each box contains one half of the balls, up to fluctuations, and the 
configuration where box number 1 is empty. This entropy difference plays the role 



of a reduced activation energy in (1^). This is an elementary example of an entropy 
barrier. 
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1.2. The Monkey urn model 

Let us introduce a variant of the Ehrenfest urn model, which we choose to call the 
Monkey urn model, because it corresponds to the image of a monkey playing at 
exchanging balls between two boxes. The key difference with the previous case is 
that now at each time step a box (either 1 or 2) is chosen at random, and one ball is 
moved from the chosen box (provided it is non-empty) to the other one. The choice 
of a box, instead of a ball, induces drastic changes in the statics and dynamics of the 
model. 

The equilibrium state of the Monkey urn model is now characterized by uniform 
occupation probabilities 

1 

N + 1 

The evolution of the occupation probabilities is governed by the master equation 

describing a usual random walk over the integers 0, . . . , A^, with appropriate reflecting 
boundary conditions. 

The relaxation rates of the model read Am = 2 sin^ (TO7r/(2(A^ + 1))), with 
m = 0, . . . , N. The relaxation time reads therefore 

Teq^^^^iV^, (1.6) 
Al TT^ 

where the scaling r^q ^ N"^ originates in the diffusive nature of the motion of balls. 



/fe,oq = ^cq(iVl = fc) = (fc-0,...,iV). (1.5) 



Because of the flatness of the equilibrium distribution (1.5), there are no entropy 
barriers in the Monkey urn model. In particular, the mean time Tq for a box to get 
empty scales as iV^, proportionally to the relaxation time (1.6). 



2. Many-urn models: statics and dynamics 



We now turn to the discussion of static and dynamical properties of two classes of 
many- urn models, generalizing the Ehrenfest and the Monkey urn models described 
so far 1^. This classification encompasses examples of urn models studied in recent 
years. The backgammon model is a prototype of the Ehrenfest class, while model 
B of References [|[ ^ and the zeta urn model ^ belong to the Monkey class. 

A many-urn model consists of N balls, distributed among M boxes. Box number i 
contains Ni balls, with 

M 

1=1 

Most investigations of many-urn models have dealt with the thermodynamic limit 
{M — > oo, — > oo), at a fixed density p ~ N/M of balls per box. 



2.1. Statics 

Two ingredients are necessary in order to define static properties of many-urn models 
within the framework of statistical mechanics: 
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A priori statistics. To each configuration {N^} is attributed an a priori statistical 
weight 

f 1 

Wm)^l]^^m (Ehrenfest class), ^^^^^ 

I 1 (Monkey class). 

The formula for the Ehrenfest class generalizes the binomial law (|1 . 1|) , and corresponds 
to the Maxwell-Boltzmann prescription for indistinguishable particles. The formula 
for the Monkey class generalizes the flat distribution (1.5), and formally corresponds 
to Bose- Einstein statistics for the occupation numbers of quantum states . 

Hamiltonian. The Hamiltonian of a many-urn model is defined as the sum of inde- 
pendent contributions of boxes, of the form 

M 

^({^^}) =E^(^»)' (2.2) 

i=l 

so that the unnormalized Boltzmann weight of a configuration {Ni} at temperature 
T= 1/13 is 

M 



with 



e ^Wpn^. (2.3) 



The canonical partition function thus reads 

Z{M,N) = e-^''"'^*" j(^iY„7v] . 



Using (^j|) and ( |2.3| ), and a contour-integral representation of the Kronecker symbol, 
we are left with 

m,N) = j,^r^^[P{z)t\ (2.5) 

where 

P(7\ - V j-fc X / (Ehrenfest class), , . 

^^^^ - Z^^fe^ ^ I 1 (Monkey class). ^^■^> 

Static properties of many-urn models are therefore entirely encoded in the 
temperature-dependent generating series P{z) of the Boltzmann weights pk of (2.4). 
In the thermodynamic limit, at fixed density p, the free energy per box, 

I3F^- lim -i- In Z{M, N) {N w Mp), 



can be obtained by evaluating the contour integral in ( |2.5| ) by the saddle-point method. 
The saddle-point value Zg depends on temperature and density through the equation 

'.£lM-p (27) 
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and the free energy per box reads 

PF = plnZs-lnPiz,). (2.8) 

The equilibrium occupation probabilities fk,cq of box number i do not depend 
on the box under consideration. (In the sequel we will refer to box number 1 for 
definiteness.) We have 

_ PkZ^ Jl/fc! (Ehrenfest class), , . 

•''^'^^ " { 1 (Monkey class). ^ 

The above results further simplify at infinite temperature, where pk — 1, 
independently of the form of the energy function. Equations (|2.5|)-(|2j|) become 



P(z) = c^ Z(M,iV)-— , zs^p, 

k 

l3F^p\np-p, fk^,^^e-P^ 



(2.10) 



fc! 



for the Ehrenfest class, and 



^ ' 1- z' ^ ' ^ 7V!(A-f - 1)! 



(2.11) 



/3F-plnp-(p+l)ln(p+l), A,eq ^ (^+1)^+^ 



for the Monkey class. 



2.2. Dynamics 

An elementary step of the dynamics consists in moving a ball from a departure box d to 
an arrival box a. A move is proposed to the system at every time step St = 1/M. The 
precise definition of a dynamical many-urn model should state the following points: 



Geometry. The departure box d and the arrival box a have to be connected to each 
other. Most investigations of dynamical urn models have dealt with the mean-field 
geometry, where any two boxes are connected to each other. 

Dynamical urn models can also be defined on regular lattices in finite dimension, 
with the convention that urns are located at the sites of the lattice, and that only 
nearest neighbors are connected. 



A priori dynamical rule. The proposed moves are chosen according to some a priori 
dynamical rule. There are three natural rules 

• Rule A (ball-box): Choose a ball (this defines the departure box d). Choose an 
arrival box a, connected to d. 

• Rule B (box-box): Choose a departure box d. Choose uniformly a ball in d. 
Choose an arrival box a, connected to d. 

• Rule C (ball-ball): Choose a first ball (this defines the departure box d). Choose 
a second ball (this defines the arrival box a). Propose the move if the boxes d 
and a are connected. 

The urn models described in Section 1 suggest what should be the right relation- 
ship between a priori statistics and dynamical rule. It indeed turns out that Rule A 
corresponds to the Ehrenfest class, while Rule B corresponds to the Monkey class. 
Rule C corresponds to yet another class of dynamical urn models, with no obvious 
static counterpart 
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Thermal rule. The following thermal rules respect detailed balance: 

• Metropolis: A proposed move is performed with probability 

min(l,exp(-/3A7^)), (2.12) 

where AH = H({A^/}) — H{{Ni}) is the energy difference between the configura- 
tions {Ni} before and {N^} after the proposed move. 

• Heat-bath: Among all the possible a priori choices for the arrival box a, one 
is chosen with an appropriately normalized probability, proportional to the 
Boltzmann weight exp ( — /9H({iV/})) of the configuration after the move. 



Throughout the following, we restrict ourselves to the mean-field geometry, and 
to a homogeneous disordered initial state, such as, e.g., the infinite-temperature 
equihbrium state ( 2.10|) , ( ^.ll| ). 

The key quantities are again the occupation probabihties fk{t) — V{Ni{t) = k). 
Their evolution is given by the master equation 

= ^(^fc+i,Ki)+7r£,fc-i(t)-^fc,K0-^^,fc(i)), (2.13) 

where nk,e{t) denotes the contribution of a move from the departure box d, containing 
Nd = k balls, to the arrival box a, containing Na = £ balls. 
The a priori dynamical rule yields 

k (Ehrenfest class), 
1 (Monkey class). 

The factor (1 — 6k,o) accounts for the fact that an actual move only takes place if the 
departure box d is not empty. The acceptance rates Wk/(t) depend on the thermal 
part of the rule. With the notation (2^), the rates for the MetropoHs rule [^, 

Wk,e = mm 1 



TTkAt) = fkit)fi{t)WkAt)i^ - 4,o) X 



PkPe+i , 

are independent of time, while those for the heat-bath rule 



only depend on the la bel £ o f the arrival box. 

In analogy with (|l|), ( |2.13| ) can be recast as the master equation of a random 
walk over the positive integers (fc = 0, 1, . . .): 

- Mfe+i W/fe+i W + A,-i(OA-i(t) 

- (A*fe(i) + Afc(t))/fe(t) (fc>l), (2.14) 
"^^"^^^ -/.i(t)/i(t)-Ao(i)/o(t), 



At 



where 



Afc(0 = -7^E'^"W' A^feW^TTirE'^MW (2-15) 
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are, respectively, the hopping rate to the right, corresponding to a = 1, i.e., Ni — k 



Ni — k ^ 1, and to the left, corresponding to c? = 1, i.e., A^i 
equation for /o(i) is special, as A_i(t) = /xo(t) = 0. 

The master equation preserves the sum rules 



k>Q 



A(t) = i, J2kMt)^{N,{t))^ 



1. The 



(2.16) 



fe>0 



expressing the conservation of probability and of the number of balls. 

The rates Afc(t) and /i.fe(t), defined in (2.15), depend on the occupation probabili- 
ties fk{t) themselves, so that the master equation (2.14) is actually a nonlinear set of 
first-order equations for the fk{t). The equilibrium occupation probabilities ( |2.9| ) are 
recovered as their unique stationary state. 



3. Equilibrium vs. nonequilibrium dynamics: two-time quantities 

Nonequilibrium properties of dynamical urn models such as the backgammon model |0, 
H, |ll|, |l^ or the zeta urn model ^ have been extensively studied in recent years. 

These models are simple enough to allow analytical investigations, and rich 
enough to illustrate most features of nonequilibrium dynamics referred to as aging. 
The latter include the scaling properties of the correlation and response functions 
of various observables, and of their fluctuation-dissipation ratios. In this section we 
recall how the temporal evolution of these two-time quantities can be derived within 
the master-equation formalism of Section 2. 

We consider a dynamical urn model of either class, in the thermodynamic limit. 
For definiteness, our observable will be the local density of balls, defined as the 



fluctuating population Ni{t) of box number 1, say. The second of the sum rules (2.16) 
ensures that the mean density (iVi(i)) = pis constant in time. Its fluctuations Ni[t)~p 
exhibits, however, non-trivial dynamical properties, that can be characterized by two- 
time correlation and response functions, to be defined below. Hereafter, time variables 
will be such that Q < s < t. The earlier time s is referred to as the waiting time, the 
later time t as the observation time. 

3.1. Correlation function 

The (local, connected) two-time correlation function of the density is defined as 

C{t,s) = {N,{s)N^{t))^p\ 

This can be recast as 

C(t,s) = ^fc7fe(t, s)-p', 
fe>i 

where the quantities 

lk{t,s) = ^j7,(5)P{A^i(i) = k I N,{s)^j} 

assume the initial values 



7fe(s, s) = kfk{s) 
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at t — s. Their temporal evolution for t > s is given by a master equation similar 

= /i/c+i {t)"fk+i [t, s) + Afc-1 {t)"fk-i [t, s) 



to ( |2.14D : 



dt 

-{lik{t) + \k{t))ik{t,s) (fc>l), (3.1) 

g^' ''^'^ = A*i (^71 {t,s)~ Ao (i)7o (i , s) . 

The rates \u{t) and /^^(i), defined in ( p.l5| ), only depend on the /^(i), so that (3.1) 
are linear equations for the 7fc(i, s). 

3. 2. Response function 

The (local) response function measures the influence on the mean density in box 
number 1 of a perturbation in the canonically conjugate variable, i.e., the local 
chemical potential acting on the same box. 

Suppose that box number 1 is subjected to a small time-dependent chemical 
potential cki(i), so that the total Hamiltonian of the system is now 

M 

^({^0) =5I^(^^)+"lW^l■ 
The mean density in box number 1 reads 

t 



{Ni{t)) = p+ R{t,s)ai{s)ds- 





where only the term linear in a(s) has been written explicitly, and causality has been 
used. The kernel of the linear response: 

S{N,{t)) 

is called the two-time response function. The temporal evolution of this function can 
again be determined by means of the master-equation formalism. 

3.3. Fluctuation-dissipation theorem: equilibrium vs. nonequilibrium behavior 

Dynamical urn models, as any statistical-mechanical model, have a fast convergence to 
equilibrium, with a finite relaxation time Teq, in their high-temperature (disordered) 
phase. If the earlier time exceeds the relaxation time (s 3> T^q), the system is at equil- 
ibrium. One-time quantities assume their equilibrium values. Two-time quantities, 
such as the correlation and response functions, are invariant under time translations: 

C{t, S) = Ccq(r), R{t, S) = i?eq(T), (3.2) 

where 

T = i - s > 0. 

Furthermore, these quantities are related by the fluctuation-dissipation theorem 

i?c,(r)^-i^. (3.3) 

Some many-urn models exhibit a more interesting nonequilibrium behavior at 
lower temperatures, where the relaxation time Teq becomes either very large or infinite, 
so that the waiting time s and the observation time t can be much smaller than Teq. 
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In such a nonequilibrium situation, both time-translation invariance (3.2) and 
the fluctuation-dissipation theorem (3.3) are violated in general. It has become 
customary [ p^ p^ , [T^ , to characterize departure from equilibrium by the fluc- 
tuation-dissipation ratio 

R{t,s) 



Xit, 



= T 



dCjt, 
~~d7 



(3.4) 



In general, this dimensionless quantity depends on both times s and t and on the 
observable under consideration. It often exhibits a non-trivial scaling behavior in the 
two-time plane, which is viewed as one of the salient features of aging. 

The definition (3^) can be rephrased in terms of an effective temperature fl^ : 

1 dC{t,s) _ . T 



R{t, 



Tcffit, s) 



Teff(t,s) ds ' ' X{t,s)' 

At equilibrium, the fluctuation-dissipation theorem is valid, so that X{t,s) — 1 and 
Tgff (t, s) = T. For a nonequilibrium situation started from a disordered initial state, 
such as (2.10), (2.11), the effective temperature is expected to lie somewhere between 
the formal initial temperature T(0) = oo and the final one, T{oo) = T, i.e., 

< X{t,s) < 1. 

These inequalities are indeed obeyed in all the known cases. 



4. The backgammon model 

The backgammon model is a simple example of a system which exhibits slow relaxation 
due to entropy barriers |^] . This mean-field many- urn model belongs to the Ehrenfest 
class, with Metropolis dynamics, and the following choice of an energy function: 

E{N^) = -J(iV„0). 

The Hamiltonian ( ^.2| ) is therefore equal to minus the number of empty boxes, so that 
balls tend to condensate in fewer and fewer boxes as times passes, at least at low 
temperature. 



Statics 

The equilibrium occupation probabilities read 



The parameter Acq, equal to in the formalism of Section 2, represents the mean 
population of a non-empty box. It is related to temperature and density by 



implying 



Aeq«/3-ln^ 
P 



at low temperature. The roughly linear divergence of Aeq is the signature of a zero- 
temperature condensation transition. The density plays no particular role in this 
phase transition, so that we set p = 1 in the sequel. 
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The slow energy relaxation at low temperature can be understood, at least 
qualitatively, in terms of the concept of entropy barrier, introduced for the Ehrenfest 
model. The late stages of the dynamics indeed consist in emptying boxes whose typical 
population is Acq. The formula (1.3) suggests the scaling law InTeq ~ Aeq ~ (3. The 
relaxation time of the energy is indeed known |12[ to read 



' cq 



^/(A 

Acq 



cq/ 



e 
A, 



cq 



1 



1 



A, 



2 In /3 + 1 



(4.1) 



at low temperature, with the other characteristic times remaining of order Acq. 
function 



The 



/(A) 



dz- 



1 



Ei(A) -lnA-7, 



/o z 

where Ei is the exponential integral and 7 is Euler's constant, plays a central role in 
the analysis of the low-temperature regime. 



4-2. Nonequilibrium dynamics: logarithmic relaxation 

In the present case, the master equation ( ^.14| ) for the occupation probabilities reads 
dfkit) 



dt 

dm 

dt 

dfojt) 
dt 



k + 1 
A(^ 
2 



fk+iit) + fk-iit) 



fkit) (fc>2), 



f2{t)+^l{t)foit)~2h{t), 



A(t) 

^fiit)-l^{t)fo{t) 



(4.2) 



with 



= 1 - (1 - C-'3)/o(i), fiit)=C-^ + il 



The parameter A{t) again represents the mean population of a non-empty box in 
the low-temperature regime, while the energy density reads 

Eit)^-m^-i + ^^. 

The nonequilibrium regime of interest, called the alpha regime in the context of 
the physics of glasses, corresponds to low temperature and long times. It is natural 
to investigate the slow relaxation dynamics in this regime by means of an adiabatic 
decoupling, similar to the Born-Oppenheimer scheme, between the slow evolution of 
A(t), due to emptying new boxes, generalizing the time scale r^q, and the fast evolution 
of the other degrees of freedom, corresponding to the diffusive motion of balls between 
non-empty boxes, generalizing the relaxation times of order Acq. 

This decoupling scheme can be implemented in several ways. An analytical 
approach, extensively used in ]l^ , is based on the asymptotic analysis of formal 
integral representations of the solution to the master equation (4.2). It amounts to 
neglecting only exponentially small terms, either in A or in Acq. This approach shows 
that A{t) follows an effective Markovian dynamics 

dA 1 / _ (Ajt) ~ l)e^(*) 

dt ~ I{A{t)) V (Acq - l)eAe, 



(4.3) 
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throughout the alpha regime {A{t) and Aeq both large) 

arei 

In 

A{t) w lni + lnlnt + 



At zero temperature, the large parenthesis in (4.3) equals unity, so that 

Inlnt - 2 



Int 

The same logarithmic law is present at finite temperature in the aging regime, namely 
for A(t) <C Acq, i.e., t <C Toq. It is the nonequilibrium counterpart of the equilibrium 
relationship (4.1) between r^q and Acq. 



4-3. Nonequilibrium dynamics: two-time quantities 

The two-time quantities defined in Section 3 have been studied for both local energy 
and density fluctuations, by means of the above adiabatic approach [T^ . The main 
results are as follows. 

Throughout the alpha regime, the two-time correlation and response functions 
exhibit a leading power-law scaling, with prefactors which depend logarithmically on 
both times s and t. Consider for definiteness zero-temperature dynamics. The energy 
and density correlation functions read 

^.(M)-|@, C.(M)-(A(.)-l)-^(^)*^^) 



$(t) ' ' ' ' ' 'A(s)$(i)' 



with 



$(i) = exp 

= Vt Int,. , 

int 

so that their leading long-time behavior is 




C.(t,.)«^^, C.(M)-y^ln.. 

Similar scaling expressions can be obtained for the local energy and density 
response functions Rsit^s) and R^(t,s). It turns out that two density response 
functions have to be considered, as a peculiarity of Metropolis dynamics, because the 



acceptance rate (2.12) is not differentiable (it has a cusp) at AH — 0. The associated 
fluctuation-dissipation ratios are independent of the observation time t throughout 
the alpha regime, provided r = i — s3>lisa macroscopic time. 
For a large waiting time s, one has 

XEit,s) « 1 
X+{t,s)^l 

^^^^''^""^ Aisr Aisf Aisr ' 
The fluctuation-dissipation theorem is thus asymptotically restored in the long- 
time regime, with a leading logarithmic correction in l/(lns)^, which depends both 
on the observable (energy vs. density fluctuations) and on microscopic details of the 
model (Metropolis vs. heat-bath dynamics). 



2 


4 


6 


Aisr 


Msf 


A(.)^ + 


2 


2 


12 


A(.)^ 


A{sf 


A(.)^ + 


3 


3 


15 

+ 
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5. The zeta urn model 



The static zeta urn model has been introduced |]1^ as a mean-field model of discretized 
quantum gravity. Its dynamics has been subsequently investigated |j. The zeta 
urn model belongs to the Monkey class, with mean-field heat-bath dynamics, defined 
by the energy function 

E{m) = ln(7V, + 1). (5.1) 



5.1. Statics 

The logarithmic confining potential (|]^) is such that the Boltzmann weight 

falls off as a power law, with a temperature-dependent exponent. The corresponding 
generating series has a power-law singularity at Zc — 1, of the form Psg{z) ~ z^~^, 
allowing a condensation transition for finite values of temperature and density. At low 
enough temperature (/3 > 2), Psgiz) vanishes more rapidly than linearly as z — > 1, so 
that there is a finite critical density, corresponding to Zs = Zc = 1, namely 

^ P'jl) ^ C(/? - 1) - C(/3) 

p(i) cm 

where C denotes Riemann's zeta function, which also enters (x 
the model, proposed in The static phase diagram of the model is as follows |l 



, hence the name of 



Fluid phase [p < pc). In this low- dens ity phase we have Zg < 1, so that the 
equilibrium occupation probabilities ( [2.9| ) decay exponentially. The density goes 
continuously to its critical value {p ^ Pc) as — *■ 1~. 



Criticality {p = pc). Along the critical line in the temperature-density plane, the 
equilibrium occupation probabilities 

h -i^iiri (52) 

/fexq - ^^^^ (5.2) 

obey a zeta distribution, with a temperature-dependent exponent. More generally, 
critical exponents and other universal quantities depend continuously on temperature 
along the critical line. The variance of the population of box number 1: 

Coq — ^ fk.cq — pI = ^ 7Tff]2 ' i^-'^) 

k>0 ' 

is finite for /3 > 3 (regular critical regime), while it is infinite in the strong-fluctuation 
critical regime (2 < /3 < 3). 



Condensed phase {p > Pc). A macroscopic condensate appears in this high-density 
phase. Indeed, the static configurations are such that the critical occupation proba- 
bilities ( ^.2[ ) still hold for all boxes but one, in which an extensive number of balls, of 
order N — Mpc — M{p — Pc), is condensed. 
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5.2. Nonequilibrium critical dynamics (p = p^) 

Nonequilibrium properties of the zeta urn model have been studied recently Q , both 
at criticality and in the condensed phase, pursuing earlier investigations The most 
salient results are as follows. 

Consider the situation of regular criticality (/? > 3), with a disordered initial 



state such as ( 2.11 ). For any large but finite time t, the system looks critical, i.e., 
the occupation probabilities fk (t) have essentially converged toward their equilibrium 
values ( |5.2[ ), for k <C t^^^, while for k ^ f^/^ the system still looks disordered, i.e., the 
fkit) fall off very fast. At a quantitative level, this crossover behavior is described by 
the scaling law 

fk{t) ^ fk,,^F{u), u = kt-^'^. (5.4) 
The scaling function F(u) is determined by explicitly solving the differential equation 



describing the continuum limit of the master equation (2.14). Among other conse- 
quences, the variance of the population of box number 1 is found to converge to its 
equilibrium value (|5.3|) as a power law: 



93-/3 ^-(/3-3)/2 

at, t) = {N^ - pl « C., - (^,3)r((;,+ i)/2)C(/3) - (5-5) 

The two-time local density correlation and response functions have been investi- 
gated along the same lines In the nonequilibrium scaling regime {s,t ^ 1), these 
two-time quantities are found to scale as 

C{t,s) « S-('3-3)/2 $(3.)^ 

Sl£) «S-(/3-l)/2$^(^) (5.6) 
as 

R{t,s) W S-('5-l)/2 $2(X), 

where the dimensionless time ratio 

x^tls>l 

is the temporal analogue of an aspect ratio. As a consequence, the fluctuation-dissi- 
pation ratio X{t, s) only depends on x in the scaling regime: 

X{t,s)^X{x)^T^^ 
^i{x) 

The above results illustrate general predictions on nonequilibrium critical dyna- 
mics 1 18 , 19|, pO|, ^ . The exponent of the waiting time s in the first line of ( 5.6 ) already 



appears in (5.5). It is related to the anomalous dimension of the observable under 
consideration, and would read (rf — 2 + if)/ Zc for a d-dimensional ferromagnet, with 77 
being the equilibrium correlation exponent and Zc the dynamical critical exponent. 
The scaling functions '^{x), <&i,2(2^) are universal up to an overall multiplicative 
constant, and they obey a common power-law fall-off in a;^'^/^. The latter exponent is 
not related to exponents pertaining to usual (equilibrium) critical dynamics. It would 
read — Ac/^c = 6c — d/zc for a ferromagnet, where Ac is the critical autocorrelation 
exponent and 8c is the critical initial-slip exponent |l^ . 

As a consequence, the dimensionless scaling function X{x) is universal, and it 
admits a non-trivial limit value in the regime where both time variables are well 
separated in the scaling regime: 

A'oo = lim lim X{t, s) = A'(oo). 

s — *oo t — ^00 
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Explicit expressions for the above scaling functions have been derived, using a 
spectral decomposition in Laguerre polynomials. Let us mention that the limit fluc- 
tuation-dissipation ratio, 

X^^l±l (/?>3), 

lies in an unusually high range (4/5 < < 1) for a critical system. Indeed, 

statistical-mechanical models such as ferromagnets are observed to have < X^o < 1/2 
at their critical point. The upper bound X^c — 1/2, corresponding to the mean- field 
situation , is also observed in a range of simpler models [|l3[ |l^, |l^, |l^ . 



5.3. Nonequilibrium condensation dynamics {p > Pc) 



Consider now a disordered initial condition such as ( 2.11 ) with p > pc, so that the 
system will build up a macroscopic condensate. The main results concerning this 
situation are as follows 

The growth of the condensate exhibits a scaling regime for 1 ^ t ^ IVP, 
while finite-size effects become dominant for later times. In this scaling regime, 
the occupation probabilities fk{t) have essentially converged toward their equilibrium 



values (5.2) for k <C t^^^, while each box has a small probability, of order t~^^^, to be 
part of the growing condensate, and to contain a large number of balls, of order t^/^. 
This picture can again be turned to a scaling law, which is of the form 



Contrary to the scaling function F{u) of (5.4), neither the present scaling function 



nor those entering (5.7) can be expressed in closed form, because the continuum limit 



of the master equation (2.14) in the condensed phase involves a biconfluent Heun 
differential equation, whose solutions are not known explicitly. 

The investigation of the scaling regime of condensation dynamics can be extended 
to two-time quantities |^ , which are found to scale as 

C{t,s)^{p-p,)s^/^<^{x), 

^^«(p-p.).-V2<J,^(,)^ 

i?(M)«(p-p,)s-l/2ci,2(x), (^-^^ 

X(M)-A'(x)=T|iM. 

This condensation dynamics exhibits many peculiarities with respect to usual 
phase-ordering dynamics which takes place, e.g., if a ferromagnet is quenched 
below its critical temperature. In the latter case, domain growth and phase separation 
take place in a statistically homogeneous way, at least for an infinite system. To the 
contrary, in the present situation, condensation takes place in a very inhomogeneous 
fashion, as fewer and fewer boxes are involved. 

Among unusual features of two-time quantities, let us emphasize that the scaling 
functions ^{x), $1,2(3;) assume finite values, both at coinciding times {x = 1) and 
in the limit of large time separations (x = 00). Furthermore, the limit fluctuation- 
dissipation ratio Xoo — A" (00) depends continuously on temperature throughout the 
condensed phase {(3 > 2, i.e., T < 1/2), and it vanishes only as 

7.3/4 

= - i_ + . . . 
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at low temperature, while coarsening systems are known [ p^ to have identically 
= throughout their low-temperature phase. 
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